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Spacecraft  periodic-disturbance  rejection  using  a  realistic  spacecraft  hardware  simulator  and  its  associated 
models  is  investigated.  The  effectiveness  of  the  dipole-type  disturbanee  rejection  filter  on  the  current  realistic 
nonlinear  rigid-body  spacecraft  model  is  validated.  However,  it  is  shown  that  the  rejection  filter  is  not  robust 
to  disturbance  frequency  uncertainty.  Therefore,  system  identification  techniques  are  needed  if  the  disturbance 
frequency  is  unknown  or  changing.  Based  on  experience  with  the  current  test  bed,  two  simple  closed-loop  system 
identification  methods  are  introduced  to  experimentally  identify  the  disturbance  frequency  using  data  from  the 
control  input  torque  signal.  It  is  shown  that  an  incorrect  filter  frequency  results  in  beating  of  the  control  input 
torque.  Examination  of  the  beat  envelope  leads  to  tuning  of  the  rejection  filter. 


I.  Introduction 

A  COMMON  challenge  in  high-performance  spacecraft  control 
is  the  mitigation  of  periodic  disturbances  caused  by  internal  or 
external  components.  Examples  include  disturbances  caused  by  con¬ 
trol  moment  gyroscopes  (CMGs)*  and  disturbances  caused  by  the 
excitation  of  flexible  appendages,  as  was  the  case  with  the  Hubble 
Space  Telescope.^  In  the  present  paper,  disturbance-rejection  filters 
are  studied  using  the  Naval  Postgraduate  School  (NPS)  Three-Axis 
Simulator  (TAS)  test  bed.  The  TAS  was  developed  at  the  NPS  in 
Monterey,  California  and  is  shown  in  Fig.  1.  The  purpose  of  the 
TAS  was  to  prove  a  satellite’s  ability  to  redirect  a  laser  originating 
from  Earth,  an  aircraft,  or  another  satellite  to  another  location  on 
Earth  for  communication  or  defense  applications.  A  satellite  utiliz¬ 
ing  lasers  requires  exceptional  acquisition,  pointing,  and  tracking 
capabilities.  On-going  research  has  been  conducted  to  study  fine 
acquisition,  pointing,  and  tracking  requirements.^  '*  Periodic  distur¬ 
bances  can  be  introduced  into  the  test  bed  by  deflecting  a  flexible 
beam  on  the  test  bed. 

Disturbance  rejection  of  periodic  disturbances  is  a  classical  field 
that  has  been  examined  using  a  variety  of  methods,  including  both 
feedforward  and  feedback  methods.  The  current  paper  further  exam¬ 
ines  the  dipole  disturbance-rejection  filter  (DRF),  which  explicitly 
employs  the  periodic-disturbance  frequency  into  a  filter  embedded 
in  the  feedback  control  loop.  The  DRF  is  based  on  the  internal  model 
principle,^"^  and  its  attenuation  effectiveness  has  been  shown  in  a 
simulation  of  the  Hubble  Space  Telescope*  and  in  simulations  for  a 
spacecraft  with  CMGs.*'®  The  filter’s  effectiveness  was  also  shown 
experimentally  on  NASA’s  minimast  truss  structure  at  the  NASA 
Langley  Research  Center'**  and  on  NASA’s  Advanced  Control  Eval¬ 
uation  for  Structures  test  bed  at  the  NASA  Marshall  Space  Flight 
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Center."  The  Hubble  Space  Telescope  solar  array  disturbance  prob¬ 
lem  was  also  attenuated  by  the  same  filter,  although  the  filter  was 
not  termed  the  dipole  DRF.^ 

The  dipole  DRF  is  based  on  linear  theory.  Therefore,  it  is  of 
interest  to  examine  its  efficacy  on  the  current  nonlinear  spacecraft 
model.  In  addition,  the  robustness  of  the  rejection  filter  is  explored  in 
the  current  paper.  It  is  shown  that  the  proposed  rejection  filter  is  not 
robust  to  frequency  uncertainty.  Based  on  experience  with  the  TAS, 
simple  closed-loop  system  identification  methods  are  introduced  to 
identify  experimentally  the  disturbance  frequency  to  be  used  within 
the  filter  design. 

The  paper  is  organized  as  follows.  In  Sec.  II,  the  theory  behind  the 
DRF  is  reviewed.  In  Sec.  Ill,  the  hardware  test  bed  and  associated 
linear  and  nonlinear  system  models  are  described.  In  Sec.  IV,  the 
DRF  is  demonstrated  on  the  test  bed  models.  In  Sec.  V,  robustness  of 
the  filter  to  frequency  uncertainty  is  investigated.  As  expected,  the 
DRF’s  attenuation  effectiveness  is  based  on  accurate  knowledge  of 
the  disturbance  frequency.  In  Sec.  VI,  motivated  by  test  bed  results, 
two  methods  are  introduced  to  determine  experimentally  the  dis¬ 
turbance  frequency  when  using  DRFs,  and  the  advantages  of  each 
are  discussed.  Finally,  in  Sec.  VII,  the  paper  is  concluded  with  a 
summary  of  results  and  future  research  directions. 

II.  Dipole  DRF 

The  dipole  DRF  is  based  on  the  internal  model  principle,  which 
places  periodic-disturbance  poles  inside  the  compensator  loop  so 
that  the  unstable  disturbance  poles  are  canceled  during  loop  closure 
between  the  output  signal  and  input  disturbance  (Fig.  2).^  Equa¬ 
tion  (1)  shows  the  closed-loop  transfer  function  between  the  in¬ 
put  disturbance  signal  GjjistC'S)  and  output  signal  Y (v)  of  Fig.  2.  In 
closed  loop,  the  disturbance  poles  are  canceled  by  the  inclusion  of 
a  DRF, 


T(v)  = 


_ l/DC) _ 

1  -F  N(s)NAs)ID{s)DAs)Dii,ds) 


Ndhl  , 

Although  the  traditional  DRF  shown  in  Fig.  2  does  not  include 
a  numerator,  a  numerator  with  the  same  order  as  the  denominator 
should  be  chosen  to  ensure  a  zero  for  every  pole  in  the  DRF  root 
locus.®’’  To  derive  the  DRF,  first  consider  the  Laplace  transform 
of  a  persistent  periodic  disturbance  of  amplitude  A  and  a  single 
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Fig.  1  NFS’s  Three-Axis  Simulator  test  bed. 
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Fig.  3  DRF  Bode  diagram. 
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Fig.  2  Closed-loop  disturbance  rejection  control  system  based  on  in¬ 
ternal  model  principle. 


frequency  to  that  can  be  represented  as 


fr  dist  (“^  ) 


A /to 
jop-  -f  1 


(2) 


Only  the  disturbance  denominator  (pole)  is  of  interest.  Therefore, 
in  the  filter  design,  the  persistent  periodic-disturbance  frequency  to 
is  employed  as  the  disturbance  rejection  frequency  to^,  and  the  DRF 
becomes 


GDRF(i')  = 


+  1 

i^/toj  -f  1 


(3) 


Inclinometer^ 
(Pitch  Position) 
Note:  Not  Sho^ 


Be^ctiorrWheels 


■~'iilnclinometer1 
(Roll  Position) 


Rate  Gyros 
On-Board  Computer 


frequency,® 


Fig.  4  TAS  components. 
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III.  NFS’s  TAS  and  Its  Linear  and  Nonlinear 
Simulation  Models 


Note  that  the  DRF  numerator  is  included  as  a  modification  of 
Eq.  (2),  and  recall  that  the  numerator  and  denominator  should  be  of 
the  same  order.  The  DRF  Bode  plot  is  shown  in  Fig.  3,  and  some 
guidelines  in  selecting  a>^  are  provided  in  the  next  paragraph. 

An  cOp  and  pair  is  termed  a  dipole,  and  the  difference  between 
the  (Op  and  co^  values  is  termed  the  dipole  strength.^  The  settling 
time  of  the  transient  response  increases  as  the  difference  between 
the  dipole  terms  increases.  The  co^  can  be  chosen  to  be  less  than  or 
greater  than  cOp  as  shown  in  Fig.  3.  A  magnitude  magnification  in  the 
high  frequencies  results  from  co^xop,  and  a  magnitude  attenuation 
in  the  high  frequencies  results  from  co^  <cOp.  The  net  magnitude  at¬ 
tenuation  or  magnification  increases  as  the  dipole  strength  increases. 
It  is  suggested  that  co^  be  selected  to  be  between  two  consecutive 
poles  of  the  system  for  stability.® 

Shown  in  Eq.  (3),  the  DRE  is  basically  a  second-order  generalized 
filter^  Gcis)  [Eq.  (4)],  without  damping  terms.  Actually,  damping 
terms  can  be  included  as  reported  in  the  literature.*  In  a  Hubble 
Space  Telescope  solar  array  disturbance  analysis,  the  numerator 
damping  term  was  made  high  compared  to  the  denominator  damp¬ 
ing  term,  which  made  the  antiresonance  practically  nonexistent.* 
The  filter  was  designed  by  adjusting  the  damping  terms  to  create  a 
second-order  non-minimum  phase  bandpass  filter.  Note  that  for  dis¬ 
turbances  of  several  frequencies,  dipoles  may  be  included  for  each 


The  DRE  and  the  underlying  internal  model  principle  are  all  based 
on  linear  system  theory.  Studies  of  DRE  filters  on  linear  system 
models  of  spacecraft  have  been  reported,®’^’'^  and  results  validate 
the  linear  system  theory  as  expected.  In  this  paper,  the  DRE  is  only 
demonstrated  on  a  nonlinear  model  of  the  experimental  test  bed. 
However,  the  proposed  input  disturbance  frequency  determination 
methods  discussed  in  Sec.  VI  are  analyzed  using  both  a  linear  and 
nonlinear  simulation  model. 

The  TAS  hardware  simulator’s  components  are  shown  in  Eig.  4. 
The  simulator  underside  has  a  semispherical  ball  that  mates  with 
an  air  bearing’s  semispherical  cup.  The  TAS  essentially  has  three- 
axis  frictionless  motion  when  pressure  is  supplied  through  the  air 
bearing,  which  produces  a  thin  film  of  air  between  the  semispherical 
ball  and  cup. 

Removable  hanging  masses  and  fine-tuning  masses  are  utilized 
to  balance  the  TAS  such  that  the  center  of  mass  corresponds  to  the 
TAS  center  of  rotation.  In  doing  so,  the  TAS  emulates  a  satellite  in 
a  space  environment.  However,  perfect  balancing  cannot  be  realisti¬ 
cally  achieved  by  manual  tuning,  and  so  a  control  input  torque  must 
continually  counter  the  imbalance  disturbance  torque. 

Eor  actuators,  the  TAS  has  three  reaction  wheels  that  are  orthog¬ 
onal  to  each  other.  The  TAS  uses  two  inclinometers  and  a  one- 
axis  infrared  (IR)  sensor  for  attitude  sensing.  The  inclinometers  are 
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Fig.  5  Closed-loop  block  diagram  of  rigid-body  spacecraft  with  PID  controller,  ROF,  and  DRF. 


capacitive  liquid-based  sensors  where  the  sensors’  capacitance  is  a 
linear  function  of  the  liquid  angle.  Although  inclinometers  are  not 
used  on  actual  spacecraft,  inclinometers  provide  a  reasonably  ac¬ 
curate  and  inexpensive  alternative  to  determine  the  TAS  pitch  and 
roll  angles.  The  IR  sensor  has  two  phototransistors  that  measure 
the  change  in  intensity  of  a  light  source  to  determine  the  TAS  yaw 
angle.  For  angular  rate,  the  TAS  uses  three  orthogonal  mechanical 
rate  gyros.  The  TAS  uses  MATLAB/Simulink™  and  xPC  Target™ 
for  real-time  control.  More  details  of  the  hardware  design  may  be 
found  in  Ref.  12. 

The  TAS  is  physically  subjected  to  a  periodic  disturbance  pro¬ 
duced  from  the  initial  displacement  of  an  aluminum  point  mass 
and  beam  that  is  attached  to  the  edge  of  the  TAS  as  shown  in 
Fig.  1.  The  disturbance  parameters  were  experimentally  derived'^ 
as  A  =  2.13  N-  m  (initial  disturbance  amplitude),  cl)„  —  0.61512  Hz 
(natural  frequency),  and  f  =  0.002303  (damping  ratio).  Note  that 
in  the  current  study,  the  beam  is  considered  only  as  a  source  of 
disturbance  to  rigid-body  models. 
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Figure  6  shows  the  three-axis  nonlinear  rigid-body  TAS  Simulink 
model  used  to  analyze  the  dipole  filter.  An  attitude  command  is 
prescribed  that  is  then  transformed  to  a  rate  command  using  Eq.  (9), 
assuming  a  2-3-1  rotation  sequence.^ 
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A.  TAS  Single-Axis  Linear  Simulation  Model 

The  single-axis  rigid-body  linear  TAS  simulation  block  diagram 
is  shown  in  Fig.  5 .  The  simulation  consists  of  a  proportional  integral 
derivative  (PID)  controller,  roll-off  filter  (ROF),  DRF,  disturbance 
signal,  and  a  rigid-body  spacecraft  plant.  The  transfer  function  for 
the  PID  controller,  ROF,  and  plant  are  shown  in  Eqs.  (5-7).  The 
disturbance  signal  is  represented  in  Eq.  (2).  Similarly,  the  DRE  is 
represented  in  Eq.  (3), 

GpId(.s)  =  {(Dfc/i)[5^  +  {1/Tj)s  +  (1/TjT,)]} 


'  .Command' 
"l 
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O) 


Command 

2 


,  .Command 
CO^ 


(9) 


The  measured  or  actual  rate  is  solved  by  Euler’s  rotational  equa¬ 
tion,  which  is  shown  in  Eq.  (8).  The  rate  is  then  fed  back  so  that  a 
rate  error  signal  is  calculated  by  taking  the  difference  between  the 
commanded  and  actual  rate.  At  the  same  time,  the  actual  rate  is  used 
in  Eq.  (10)  to  calculate  the  actual  attitude  in  quaternions,’  '"* 


=  {[14.4686(1.1871)/x][i-|-27r(0.005)][i-l-27r(0.006)]}  (5) 

where  k,  t,,  and  are  the  proportional,  integral,  and  derivative 
(PID)  gains,  respectively. 

GrofCv)  =  {27r(0.9)/[x  -f  27r(0.9)]}  (6) 

G,(s)  =  [1/75’]  =  [1/555’]  (7) 

where  J  (kg  •  m’)  is  the  rigid-body  inertia  about  the  single  axis. 

The  PID  controller  and  ROF  were  designed  by  classical  frequency 
response  shaping  methods  to  maximize  the  gain  and  phase  margins 
and  were  validated  on  the  hardware  test  bed.'’  The  inertia  parameter 
shown  in  Eq.  (7)  was  experimentally  determined.*’ 

B.  TAS  Three-Axis  Nonlinear  Rigid-Body  Simulation  Model 
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The  quaternion  vector  is  transformed  to  Euler  angles  where  the 
actual  attitude  is  fed  back  to  the  commanded  attitude  so  that  the 
attitude  error  signal  is  calculated.  The  conversion  from  quaternions 
to  Euler  angles  for  a  2-3-1  rotational  sequence  is  shown  in  Eqs.  (11- 
13)  (Ref.  7).  The  attitude  and  rate  error  signals  are  utilized  in  the 
PID  controller.  A  control  input  torque  (CIT)  is  generated  after  the 
DRF  (Fig.  6). 


The  three-axis  nonlinear  rigid-body  TAS  simulation  is  based  on 
Euler’s  rotational  equation  and  the  standard  kinematic  differential 
equation  for  quaternions,  which  are  shown  in  Eqs.  (8)  and  (10), 
respectively.’’ Note  that  in  our  model,  the  roll  and  pitch  axes 
(axes  1  and  3)  are  aligned  from  the  central  point  of  the  test  bed 
table  through  the  centerlines  of  the  pitch  and  roll  reaction  wheels 
shown  in  Fig.  4.  The  yaw  axis  (axis  2)  is  aligned  with  the  vertical 
air  bearing  support  structure  shown  in  Fig.  1 . 


02  =  tan  * 


2(^1  g3  -  9'2g4) 

1  -  2{qj  +  ql)  _ 


6»3  =  sin  ^[2(qiq2  +  q3qA)} 
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(12) 


01  =  tan  * 
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(13) 


—  (Tt^'*/'")  -f  X  Tcu"/"  =  "M.on.rol  +  (8) 

dt 


The  total  external  torque  shown  on  the  right-hand  side  of  Eq.  (8) 
includes  the  CIT,  TAS  imbalance  disturbance  torque,  and  a  periodic 
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Fig.  6  Three-axis  nonlinear  rigid-body  TAS  simnlation  model. 


disturbance  torque  about  the  yaw  axis.  All  of  the  torques  are  in 
the  spacecraft  body  frame  (B)  and  are  represented  in  the  Simulink 
model  shown  in  Fig.  6.  Again,  note  that  flexibility  in  the  beam  is 
not  considered  in  the  current  model. 

The  TAS  imbalance  disturbance  torques  in  Eq.  (8)  were  approxi¬ 
mated  using  the  hardware  test  bed’s  CIT  experimental  data  and  each 
are  on  the  order  of  0.1  N  •  m.  The  imbalance  torque  is  the  result  of 
the  center  of  mass  not  coinciding  with  the  center  of  rotation.  This 
effect  changes  slightly  depending  on  the  success  of  each  balancing 
of  the  test  bed.  A  bias  in  the  CIT  experimental  data  approximates  the 
magnitude  of  imbalance  being  countered,  shown  later.  For  the  inertia 
matrix  shown  in  Eq.  (8),  the  inertia  about  the  yaw  axis  was  experi¬ 
mentally  determined.'^  The  inertia  about  the  pitch  and  roll  axes  were 
analytically  determined  as  described  in  Ref.  12.  The  products  of  in¬ 
ertia  were  assumed  small  compared  to  the  diagonal  inertias.  Small 
values  were  estimated  so  that  coupling  between  the  axes  would  also 
be  small,  based  on  the  assumption  that  our  coordinate  system  is 
close  to  the  principal  axes.  To  verify  numerically  the  final  three- 
axis  simulation  (3-AS),  the  sum  of  all  torques  in  the  zero  command 
case  was  verified  as  zero  as  expected  from  theory.*^ 


IV.  Effect  of  DRF  on  TAS  Nonlinear  Model 

The  DRF  used  in  the  simulation  experiments  is  expressed  as 
shown  in  Eq.  (14).  The  selected  follows  the  guideline  of  be¬ 
ing  between  two  consecutive  poles  of  the  system.®  The  co^  is  placed 
between  one  of  the  plant  poles  at  the  origin  and  the  DRF’s  pole. 


Gdrf(.s)  = 


+ 1 
+  1 


O.lOOi^-i-  1 

0.06695^2  -p  1 


(14) 


Figures  7  and  8  show  the  DRF’s  effectiveness  on  a  persistent 
periodic  input  disturbance  on  the  yaw.  The  attitude  command  is  a 
1-deg  position  command  in  the  yaw  and  a  0-deg  position  command 
in  the  pitch  and  roll  to  emulate  the  TAS  stabilization  process.  The 
DRF  is  extremely  effective  in  reducing  the  effect  of  the  disturbance 
on  the  yaw  position. 

Note  that  we  were  not  able  to  implement  the  DRF  on  the  hard¬ 
ware  test  bed  as  a  result  of  actuator  saturation.  Our  reaction  wheel 
actuators  were  limited  to  0.2  N  •  m,  whereas  our  required  torque  was 
2.13  N  •  m  (to  duplicate  the  simulation  results).  A  0.2-N  •  m  ampli¬ 
tude  disturbance  produces  inconsequential  effects  of  the  test  bed 
yaw  position,  due  to  the  comparatively  large  inertia  of  the  test  bed. 
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Fig.  7  DRF  effect  on  persistent  periodic  disturbance  of  0.6151  Hz  and 
2.13  N  ■  m  peak  magnitude  in  3-axis  non-linear  simulation. 


V.  DRF  Robustness  to  Frequency  Uncertainty 

Figure  7  shows  the  DRF  capability  in  attenuating  periodic  input 
disturbances.  However,  the  results  are  based  on  the  assumption  that 
the  designer  knows  exactly  the  disturbance  frequency.  In  Fig.  8 
the  response  when  the  DRF  is  designed  to  reject  the  correct  and 
several  incorrect  input  disturbance  frequencies  is  compared.  Note 
that  the  DRF  is  not  very  robust  because  there  is  negligible  attenuation 
when  the  DRF  rejection  frequency  [cUp  in  Eq.  (3)]  is  incorrect  by 
±0.033  Hz  (approximately  5%  frequency  uncertainty).  In  the  next 
section,  we  will  describe  two  simple  system  identification  methods 
to  determine  explicitly  the  filter  design  frequency. 

VI.  Input  Disturbance  Frequency  Determination 

It  was  shown  in  the  3-AS  that  the  DRF  is  not  robust  to  fre¬ 
quency  uncertainty.  If  the  input  disturbance  frequency  is  uncertain  or 
changing,  the  DRF’s  attenuation  may  be  negligible.  Therefore,  the 
ability  to  determine  the  input  disturbance  frequency  is  essential  for 
the  designer.  In  this  section,  two  closed-loop  methods  to  determine 
the  input  disturbance  frequency  will  be  shown,  based  on  experience 
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Fig.  8  Effects  of  correct  and  incorrect  DRF  rejection  freqnency  where  persistent  periodic  disturbance  is  0.6151  Hz  in  3-AS  (same  conditions  as  in 
Fig.  7). 


Fig.  9  Yaw  CIT  when  DRF  rejection  frequency  is  not  equivalent  to 
disturbance  frequency  in  single-axis  linear  simulation  (persistent  dis¬ 
turbance  frequency  =  0.6151  Hz). 


with  the  TAS.  The  first  method  evaluates  the  CIT  plot  and  simply 
counts  the  number  of  cycles  in  a  time  range.  The  second  method 
exploits  a  beating  phenomenon  in  the  CIT  that  results  from  an  inac¬ 
curate  filter  design  frequency.  Both  methods  will  be  demonstrated  in 
linear  and  nonlinear  simulations,  as  well  as  in  hardware  experiments. 

Figure  9  shows  the  CIT  plot  from  the  single-axis  linear  simula¬ 
tion  (Fig.  5).  Figure  9a  shows  the  CIT  in  closed-loop  from  the  linear 
model  when  the  DRF  is  designed  using  a  0.555 1-Hz  disturbance 
rejection  frequency  [ojp  in  Eq.  (3)],  and  Fig.  9b  shows  the  CIT  in 
closed-loop  when  the  DRF  is  designed  using  a  0.535 1-Hz  distur¬ 
bance  rejection  frequency  ojp,  where  the  actual  input  disturbance 
frequency  for  both  is  0.6151  Hz.  For  Figs.  9a  and  9b,  is  0.5  Hz 
[Eq.  (3)]. 

For  the  first  frequency  determination  method,  the  input  distur¬ 
bance  frequency  is  found  by  measuring  the  high-frequency  oscil¬ 
lation  in  Fig.  9.  This  high-frequency  oscillation  in  Figs.  9a  and  9b 
both  reveal  the  input  disturbance  frequency  of  0.6151  Hz. 

Shown  in  Fig.  9,  note  that  a  beating  phenomenon  is  occurring, 
the  result  of  a  superposition  of  two  waves  with  slightly  differ¬ 


Fig.  10  Yaw  CIT  when  DRF  rejection  frequency  is  not  eqnivalent  to 
disturbance  frequency  in  3-axis  non-linear  simulation  (persistent  dis¬ 
turbance  frequency  =  0.6151  Hz). 


ing  frequencies.'^  Through  linear  and  nonlinear  simulation  stud¬ 
ies  and  hardware  test  bed  experiments,  it  was  found  that  if  the  fil¬ 
ter  design  frequencies  obey  the  relation  to  the  actual  disturbance 
frequencies  as  <  cUp  <  actual-disturbance-frequency  in  Eq.  (3), 
then  beats  appear  in  the  Eig.  9  plot.  For  Fig.  9,  recall  that  ojp  is 
0.5551  Hz  and  0.5351  Hz  for  Figs.  9a  and  9b,  respectively.  Fig¬ 
ure  9b  shows  more  beats  than  Fig.  9a  and  also  shows  a  decrease 
in  CIT.  Therefore,  as  Wp  approaches  aj^,  which  is  0.5  Hz,  more 
beats  occur  in  the  CIT  and  the  CIT  decreases.  If  cOp  reached  qj^, 
then  Eq.  (3)  would  equate  to  1.0  as  if  the  rejection  filter  were  not 
used.  As  ajp  approaches  the  actual  input  disturbance  frequency, 
less  beats  occur,  and  the  CIT  increases  to  approach  the  disturbance 
amplitude.'^ 

To  validate  the  beating  phenomenon,  disturbances  and  filters  were 
investigated  in  both  the  three-axis  nonlinear  simulation  model  and 
on  the  actual  hardware  test  bed.  The  CIT  plot  from  the  nonlin¬ 
ear  simulation  using  the  same  disturbance  and  filter  conditions 
of  Fig.  9  is  shown  in  Fig.  10.  Note  that  beating  is  once  again 
present. 
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Fig.  11  Experimental  CIT  responses  for  0.6151-Hz  decaying  periodic 
disturbance  when  DRF  rejection  frequency  is  0.5551  Hz;  note  that  pitch 
CIT  is  biased  due  to  imbalance  of  test  bed. 


The  same  phenomenon  was  also  experimentally  validated  on 
the  hardware  TAS.  The  hardware  TAS  parameters  are  equivalent 
to  the  simulation  parameters  (which  were  originally  derived  from 
the  TAS  testbed).  Hereafter  the  TAS  attempts  to  maintain  an  Eu¬ 
ler  angle  attitude  of  [0  0  0]^  while  experiencing  a  decaying  (due 
to  beam  damping)  periodic  input  disturbance.  The  disturbance  fre¬ 
quency  is  0.6151  Hz  and  the  DRF  rejection  frequency  is  0.5551  Hz. 
The  resulting  CIT  response  is  shown  in  Fig.  11.  Again,  beats  are 
present. 

A  frequency  determination  method  employing  a  DRF  and  CIT 
beats  will  now  be  shown  using  the  nonlinear  simulation  results 
(Fig.  10).  Note  that  the  method  could  have  been  demonstrated  with 
either  the  linear  simulation  (Fig.  9)  or  the  hardware  experiment  re¬ 
sults  (Fig.  11).  The  CIT  beat  frequency  fj,,  actual  input  disturbance 
frequency  /dist,  and  DRF  rejection  design  frequency  /drf,  are  known 
to  be  related  by*^ 


fh  —  /dist  “  /drf  (15) 

Note  that  Eq.  (15)  is  simply  describing  the  definition  of  a  beat 
frequency.  Therefore,  in  this  case,  the  beat  frequency  is  produced 
by  the  slightly  differing  frequencies  of  the  actual  input  disturbance 
frequency  and  the  filter  rejection  frequency.  In  Fig.  10a,  the  beat 
frequency  (hertz)  is  shown  to  be 

/h  =  [1.5/(50  -  25)][beats/s]  =  0.06 


appear,  then  cDp  is  incorrect.  As  the  designer  alters  cDp,  either  more 
or  less  beats  will  appear.  As  stated  earlier,  if  cOp  is  increased  to¬ 
ward  the  actual  disturbance  frequency,  less  beats  will  occur,  and  the 
CIT  magnitude  will  increase  to  the  input  disturbance  magnitude.'^ 
Therefore,  by  employing  the  beats  to  determine  the  input  disturbance 
frequency,  the  designer  has  two  simple  visual  trends  that  reveal  if  the 
rejection  frequency  is  correct:  1)  decrease  in  beats  and  2)  increase  in 
CIT  magnitude.  Although  not  developed  further  here,  these  trends 
could  also  be  automated  to  tune  the  filters  on-orbit,  which  would  be 
especially  helpful  if  the  disturbance  frequency  was  changing  over 
time. 


VII.  Conclusions 

This  paper  investigated  spacecraft  disturbance  rejection  for  peri¬ 
odic  disturbances  using  the  NFS  TAS.  It  successfully  demonstrated 
a  dipole  DRF,  based  on  the  internal  model  principle.  Also,  it  showed 
that  the  filter  is  not  robust  to  frequency  uncertainty.  Designing  the 
filter  with  an  error  of  5%  in  the  disturbance  frequency  essentially 
rendered  the  disturbance  filter  ineffective.  As  a  result,  two  simple 
closed-loop  system  identification  methods  were  introduced  to  iden¬ 
tify  experimentally  the  disturbance  frequency  to  be  used  within  the 
filter  design.  One  system  identification  procedure  took  advantage  of 
the  phenomenon  that  when  a  faulty  frequency  is  used  to  design  the 
filter,  beating  appears  in  the  input  control  torque.  The  beat  frequency 
is  related  to  the  actual  disturbance  frequency. 

Three  natural  extensions  of  this  work  would  make  for  valuable 
further  work.  First,  as  noted  earlier,  we  could  not  validate  our  DRF 
methods  directly  on  the  hardware  test  bed  due  to  actuator  saturation 
in  our  reaction  wheels.  In  fact,  the  needed  reaction  wheel  torque  is 
related  to  the  disturbance  magnitude,  which  saturated  our  current 
reaction  wheels.  A  study  of  the  affect  of  actuator  saturation  on  DRF 
methods  would  be  useful.  In  this  paper,  we  considered  a  flexible 
beam  only  as  a  disturbance  source  to  a  rigid  body.  In  fact,  the  flex¬ 
ibility  of  the  beam  could  be  studied  in  two  ways:  1)  Damping  in 
the  beam  can  be  accounted  for  by  damping  terms  in  the  rejection 
filter.  2)  Possible  control-structure  interactions  could  be  explored. 
Also,  we  noted  earlier  that  our  closed-loop  methods  of  frequency 
identification  could  be  used  to  create  automatic  DRF  filter  tuning 
onboard  orbiting  spacecraft.  This  could  be  especially  helpful  if  the 
disturbance  frequency  were  changing  over  time.  The  formalization 
of  an  integrated  DRF,  system  identification,  and  system  tuning  ar¬ 
chitecture  into  an  online  adaptive  system  would  be  valuable. 
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and  in  Fig.  10b  the  beat  frequency  (hertz)  is 

fhi  =  [2/(37.5  -  12.5)][beats/s]  =  0.08 

where  /drfi  is  0.5551  Hz  for  Fig.  10a  and  /drf2  is  0.5351  for 
Fig.  10b.  Substituting  /drfi  and  f^i  intoEq.  (15)  gives  /ajs,  =  /m  -b 
/drpi  =0.6151  Hz,  and  substituting  /drf2  and  /m  intoEq.  (15)  gives 
/dist  =  fhl  +  /dRF2  =  0.6151  Hz. 

Note  that  the  first  method,  evaluating  the  high-frequency  oscilla¬ 
tion  to  determine  the  input  frequency,  could  have  been  performed 
on  any  of  the  CIT  plots.  We  would  again  find  the  input  distur¬ 
bance  frequency  to  be  0.6151  Hz.  Therefore,  if  the  filter  designer 
has  an  estimate  of  the  actual  input  disturbance  frequency,  the  de¬ 
signer  can  use  a  close  DRF  rejection  frequency  to  produce  beats  in 
the  CIT  and  explicitly  determine  the  actual  disturbance  frequency. 
Notice  that  the  simple  system  identification  method  can  be  imple¬ 
mented  whether  the  plant  is  on-site  or  remote,  such  as  on  an  orbiting 
satellite. 

When  the  two  frequency  determination  methods  are  compared, 
there  are  some  advantages  to  utilizing  the  method  employing  beats. 
Using  the  proposed  beats  method,  the  designer  is  able  to  easily 
visualize  if  the  rejection  frequency  (Dp  is  correct  or  not.  If  beats 
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